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Extensive histogram Monte-Carlo simulations ol the XY 
antilerromagnet on a stacked triangular lattice reveal expo- 
nent estimates which strongly favor a scenario of mean-field 
tricritical behavior for the spin-order transition. The corre- 
sponding chiral-order transition occurs at the same tempera- 
ture but appears to be decoupled from the spin-order. These 
results are relevant to a wide class of frustrated systems with 
planar-type order and serve to resolve a long-standing contro- 
versy regarding their criticality. 
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The nature of phase transitions in frustrated systems 
which can be mapped onto spin models has been stud- 
ied extensively over the past twenty years by means of 
renormalization-group (RG) methods and Monte-Carlo 
(MC) simulations @. The tr iangular antiferromagnet 
serves as the simplest example of geometry-induced frus- 
tration, which in this case gives rise to the non-colincar 
120° spin sturcture. This type of magnetic order can 
also be described as a helically polarized spin density. 
The appropriate Landau- Ginzburg- Wilson Hamiltonian 
is usually taken as 

H = rS ■ S* + VS • VS* + Ui(S • S*) 2 + U 2 | S • S | 2 , 

where S is a complex vector and there are two fourth- 
order terms as a consequence of frustration. One reason 
for the plethora of studies of this Hamiltonian is that in 
addition to describing helical spin systems, it is also rel- 
evant to the dipole phase of superfluid 3 He, Josephson- 
junction arrays in a transverse magnetic field, as well 
as the fully frustrated bipartite lattice (see [Q and || 
for references). Although some of the earlier RG studies 
suggested a continuous transition within standard uni- 
versality classes, others found evidence for a first-order 
transition. In view of these results, the tricritical behav- 
ior suggested by the histogram MC simulations of the 
XY stacked triangular antiferromagnet (STAF) reported 
here may not be too surprising. 

Interest in the criticality of frustrated spin systems has 
recently been enhanced due to the suggestion by Kawa- 
mura of new "chiral" universality classes associated with 
the XY and Heisenberg STAF's §,§. This claim is par- 
tially supported by arguments which demonstrate that 
the symmetry of the order parameter V in the XY case 
involves a discrete two-fold chiral degeneracy as well as 
that of the two-dimensional rotation group, so that V = 



%2 x Si y • If this is indeed the relevant symmetry, and 
if the transition is neither first-order nor tricritical, then 
the universality class should be different from the stan- 
dard ones. The strongest support for the existence of new 
universality classes comes from Kawamura's MC simula- 
tions. These were of the conventional type using rather 
large lattices Lx Lx L with L=18-60 but with a possibly 
modest number of Monte Carlo steps per site (MCS), 6- 
20 runs with 2 x 10 4 MCS each. Critical exponents were 
estimated by the conventional "data collapsing" method, 
which necessitates a simultaneous estimation of the crit- 
ical temperature. In the Heisenberg case, the reported 
values are a = 0.24(8), (3 = 0.30(2), 7 = 1.17(7) and 
v = 0.59(2). These results have recently been corrobo- 
rated by three different groups using the more accurate 
finite-size scaling based on histogram MC data |§-|7|. In 
the XY case, the exponent values reported by Kawa- 
mura are a = 0.34(6), = 0.253(10), 7 = 1.13(5) and 
v = 0.54(2). Both of these sets of exponents are quite 
different from those of any standard class in 3D. However, 
they are somewhat suggestive of mean-field tricriticality 
(apart from the values for a), where a = \, (3 — 4,7 = 1 
and v = i. Note that it is only Kawamura's estimate 
for (3 in the XY case which coincides with these values 
within error. 

Unique sets of exponents were also found for the chiral- 
ity order, which Kawamura suggests occurs at the same 
critical temperature. The coincidence, or not, of the two 
critical temperatures has received much attention in cor- 
responding 2D frustrated systems Most authors 
appear to support the notion that they are the same. 

In contrast with the suggestion of new chiral univer- 
sality, the conclusion of a study of the non-linear sigma 
model in 2 + e dimensions by Azaria et al. Jj| is that 
the criticality of such frustrated systems, at least in the 
Heisenberg case, is nonuniversal. Depending on unspec- 
ified system parameters, one can have either standard 
0(4) criticality, a first-order transition, or mean-field tri- 
critical behavior. It is natural to extend these arguments 
to the XY case and speculate either a first-order or tri- 
critical transition pQj . We note that the nature of the 
chiral transition was not addressed by these authors. 

This already confusing situation was recently further 
exacerbated by the results of Zumbach's local-potcntial- 
approximation treatment of the RG ]ll| . This work em- 
phasizes the possibility of "almost second-order phase 
transitions" for frustrated systems, where there can be 
a set of effective critical exponents (also see Refs. 0] and 
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fl2t). At least in the Heisenberg case, the distinction 
between this possibility and that of a new chiral univer- 
sality class may never be satisfactorily resolved by MC 
simulations. 

The finite-size scaling of thermodynamic functions 
evaluated from histogram MC data has demonstrated 
the ability to yield highly accurate critical-exponent esti- 
mates for unfrustrated systems |l3|Jl^ ]. In addition, this 
procedure when combined with the cumulant crossing 
method gives an independent and accurate estimate of 
critical temperatures p"3| ]. This latter feature is partic- 
ularly useful for frustrated systems in view of the pos- 
sibility that spin and chiral degrees of freedom order at 
different, but nearby, temperatures. 

It is hoped that the results of MC simulations pre- 
sented here will provide convincing evidence that the XY 
STAF exhibits mean-field tricritical behavior. 

Near-neighbor antiferromagnetic exchange coupling in 
the basal plane, J± = 1, and ferromagnetic coupling 
along the c-axis, Jii — —1, were used |Tfi|| . The Metropo- 
lis MC algorithm was employed in combination with the 
histogram technique on lattices with L= 12-33 and runs 
using 1 x 10 6 MCS for the smaller lattices and 1.2 x 10 6 
MCS for the larger lattices, after discarding the initial 
2 x 10 5 - 5 x 10 5 MCS for thermalization. Averaging was 
then made over 6 (smaller L) to 17 (larger L) runs. For 
the largest lattices, this gives a respectable 20.4 x 10 6 
MCS for averaging. The advantage of performing many 
runs is that errors can be estimated (approximately) by 
taking the standard deviation. The present work rep- 
resents one of the very few reports of finite-size scaling 
of MC data which includes error bars. All histograms 
were generated at Kawamura's estimate of the critical 
temperature, T c = 1.458. 

The correlation time r for the spin order parameter 
was estimated @ to be about 620 MCS at L = 24 and 
T = 1.458. With the assumption r ~ L 2 , it can thus 
be expected that averaging was performed using roughly 
500 independent configurations (lj] in a single run for 
our largest lattice size. Although r decreases sharply 
away from T c , it remains rather large. At T = 1.440, 
for example, r was found to be approximately 200 MCS. 
This result implies that averaging was made using not 
more than about 8 independent configurations in a single 
run for Kawamura's simulations at this temperature with 
L = 60. 

Results of applying the cumulant-crossing method flS] ] 
to estimate the critical temperatures associated with 
both spin and chiral orderings are presented in Fig. 
1. The points represent the temperatures at which the 
order-parameter cumulant U m (T) at L' crosses the cu- 
mulant at-L = 12orL=15. There is considerable 
scatter in the data and care must be taken to use only 
results with L sufficiently large to be in the asymptotic 
region where a linear extrapolation is justified |l3| ]. In the 
case of the spin order, this appears to be for Ln~ x {L' ' / L) 



<; 1.5 but somewhat larger in the case of chiral order. 
As with the Heisenberg model, finite-size effects appear 
much less pronounced for the chiral degrees of freedom 
. These results suggest that the two types of order oc- 
cur at the same temperature, T c = 1.4584(6). The possi- 
bility that there are two very close but distinct ordering 
temperatures can never be ruled-out based on finite-size 
simulations. 

The possibility that the transition is weakly first order 
was also examined. No evidence for a double-peak struc- 
tures was found in the energy histograms, consistent with 
a continuous transition. In addition, the fourth-order 
energy cumulant U evaluated at T c yielded a result ex- 
trapolated to L — > oo of U* = 0.666 652(20), consistent 
with a value | expected for a continuous transition jl8) . 
A somewhat smaller value occurs in the case of a weak 
first-order transition |l5"| . We note, however, that other 
continuous transitions which occur in this model under 
the influence of an applied magnetic field have values of 
U* closer to | than found in the present case |l5|| . Finally, 
the assumption of volume-dependent scaling of various 
thermodynamic quantities did not yield a straight-line 
fit, even for the data at large L. 

Finite-size scaling results at T c for the specific heat 
C, spin order parameter M, susceptibility (as defined in 
Kawamura's work [p|,[7|Jl3|) x> an d the first logarithmic 
derivative of the order parameter V\ = d[Ln(M)]/ dK 
(where K = T^ 1 ) G3] are shown in Figs. 2-4. Exponent 
ratios were estimated by performing Ln-Ln plots and also 
by assuming a scaling dependence F — aL x for a function 
of interest (plus a constant term in the case of the specific 
heat). Except for the specific heat where errors are very 
large, the two methods gave essentially the same results 
only if the smaller lattices L = 12 and L = 15 were ex- 
cluded. In order to estimate errors due to the uncertainty 
in T c , identical scaling was also performed at T = 1.4579 
and T = 1.4590. The resulting exponent ratios, along 
with those associated with chiral order, are presented in 
Table I. All results correspond to fits performed on data 
for L = 18 — 33, except in the case of the specific heat 
where L — 15 data was also included to reduce the er- 
ror (otherwise, the result is a = 0.47(20)). The given 
errors represent the robustness of the fitting procedure 
and do not account for error bars on the figures. Results 
for the exponents v and v K estimated from the second 
logarithmic derivative Vi Jt],[14| are 0.51(1) and 0.55(1), 
respectively. 

For ease of comparison, the results of Kawamura's 
work are also included in the table. Note that in order to 
obtain best-fit exponents by the data-collapsing method, 
Kawamura used two different values of T c : 1.458 for spin 
order and 1.459 for chiral order (within the range of our 
estimate for T c ). Within errors, however, he concludes 
that the two transitions were the same. 

From the results in Table I, it can be seen that (3, 7, /3 K , 
and 7 K are the most sensitive to the choice of T c . Vari- 
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ation in the exponents due to the error bars on the fig- 
ures was about the same. With these considerations, our 
best estimate of the exponents and their associated er- 
rors at T c = 1.4584(6) are given by: a = 0.46(10),/? = 
0.24(2), 7 = 1.03(4),!/ = 0.50(1) for spin order, and 
P K = 0.38(2), 7k = 0.90(9), v K = 0.55(1) for chiral or- 
der. These results for the spin order strongly suggest 
that the transition is mean-field tricritial. 

A possible interpretation for this behavior can be found 
in our recent results of applying an in-plane magnetic 
field p5| . This work reveals that T c is a multicritical 
point, with one phase having the symmetry of the three- 
state Potts model and consequently a weak first-order 
transition. The effect of the field is to generate a term 
third-order in S in the free energy, of the form ~ mS 3 , 
where m is the q = Fourier component of the spin den- 
sity induced by the field. Since the system is frustrated in 
the triangular plane, one might expect that short-range 
order along the c-axis is already well developed at tem- 
peratures near T c , in the present case corresponding to 
the spin component m. At zero applied field, coupling 
between this short-range order and S could generate a 
third-order term which influences the critical behavior. 
For other types of systems, there may be similar cou- 
pling to other Fourier components ]lS[| . 

It is also noteworthy that even though our simulations 
on the XY model were made with about a factor of ten 
more MCS than in the Heisenberg system JtJ , larger fluc- 
tuations were observed in the present case for the spin 
order. This can be observed by comparing the cumulant- 
crossing data of Fig. 1 with the corresponding results of 
Ref. [Q . Larger fluctuations are expected if the transition 
is at or near a tricritical point. 

Our results for the chiral order are more difficult to 
interpret. The values for the exponents are not too dif- 
ferent from those of Kawamura (when the errors are ac- 
counted for), and do not correspond to any known uni- 
versality class (also see Ref. ||). The estimates for v K are 
in agreement but are significantly different from our value 
of v for the spin order. This may indicate that the chiral 
order has a distinct correlation length and that its crit- 
icality is decoupled from the spin order fqjiol]. Another 
possibility [po[ currently under investigation is that chi- 
rality is associated with one of many possible symmetry- 
breaking fields, in a manner analogous to the infinite set 
of crossover exponents previously determined for the XY 
model Q. 

In conclusion, these results of extensive histogram MC 
simulations of the stacked triangular XY antiferromagnet 
strongly support the scenario of tricritical behavior asso- 
ciated with the spin order, in agreement with Azaria et al. 
and in contrast with the proposal of new XY-chiral uni- 
versality by Kawamura. (Recent MC simulations do sup- 
port, however, the existence of a new Heisenberg-chiral 
universality class 0-0].) This resolves a long-standing 
controversy in the litterature, and is relevant to a wide 



class of frustrated systems. Further work is necessary to 
fully understand the nature of the chiral ordering transi- 
tion which appears to be decoupled from the spin order. 
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TABLE I. Varation of exponents with assumed critical 
temperature, along with Kawamura's results. 





1.4579 


1.4584 


1.4590 


Kawamura 


a 


0.47(8) 


0.46(10) 


0.39(13) 


0.34(6) 


P 


0.23(1) 


0.24(1) 


0.26(1) 


0.25(1) 


7 


1.07(2) 


1.03(2) 


0.99(2) 


1.13(5) 


V 


0.51(1) 


0.50(1) 


0.50(1) 


0.54(2) 


A. 


0.36(1) 


0.38(1) 


0.40(1) 


0.45(2) 


7« 


0.96(2) 


0.90(2) 


0.81(2) 


0.77(5) 


v K 


0.56(1) 


0.55(1) 


0.54(1) 


0.55(2) 



FIG. 1. Results of applying the cumulant-crossing method 
(see text) to estimate the critical temperatures associated 
with the spin and chiral orderings where b — L' /L. 

FIG. 2. Finite-size scaling of the specific heat data for 
L = 12 — 33. Data at L=12 is excluded from the fit. Er- 
ror bars are estimated from the standard deviation found in 
the MC runs. 

FIG. 3. Finite-size scaling of the order parameter. Data 
at L=12 and 15 are excluded from the fit. Error bars are 
estimated from the standard deviation found in the MC runs. 

FIG. 4. Finite-size scaling of the susceptibility x as well 
as the logarithmic derivative of the order parameter Vi (see 
text), as in Fig. 2. 
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